%run preparation.py
%matplotlib inline

Dx_tip = np.diff(bmd01['x']); Dy_tip = np.diff(bmd01['y'])
DD=np.sqrt((Dx_tip**2)+(Dy_tip**2))
v=DD/(2*dt); Av_vel = np.mean(v)
vSD=np.sum(((v-Av_vel)**2)/(np.size(v)-1)); vSD=np.sqrt(vSD)
time = np.linspace(0,300,150)

fig, ax1 = plt.subplots(dpi=500)
#plt.figure(figsize=(14,10))
ax1.set_xlabel('Time', fontdict=font)
ax1.set_ylabel('Number of binding motors', fontdict=font)
ax1.plot(time,bm01, label='Change in binding motor')
#ax1.axhline(float(bm01.mean()), linestyle='--')
ax1.legend(loc='upper left')

ax2 = ax1.twinx()
ax2.set_xlabel('Time', fontdict=font)
ax2.set_ylabel('Velocity ($\mu m/sec$)', fontdict=font)
ax2.plot(time,v, 'g', label='Change in velocity')
ax2.axhline(np.mean(v), linestyle='--', color='g', label='Mean velocity')
ax2.legend(loc='upper right')
ax2.tick_params('y', colors='g')

plt.title('Binding motors [R = 0.1]', fontdict=font)
plt.grid(); plt.show()

Dx_tip = np.diff(bmd02['x']); Dy_tip = np.diff(bmd02['y'])
DD=np.sqrt((Dx_tip**2)+(Dy_tip**2))
v=DD/(2*dt); Av_vel = np.mean(v)
vSD=np.sum(((v-Av_vel)**2)/(np.size(v)-1)); vSD=np.sqrt(vSD)
time = np.linspace(0,300,150)

fig, ax1 = plt.subplots(dpi=500)
#plt.figure(figsize=(14,10))
ax1.set_xlabel('Time', fontdict=font)
ax1.set_ylabel('Number of binding motors', fontdict=font)
ax1.plot(time,bm02, label='Change in binding motor')
#ax1.axhline(float(bm01.mean()), linestyle='--')
ax1.legend(loc='upper left')

ax2 = ax1.twinx()
ax2.set_xlabel('Time', fontdict=font)
ax2.set_ylabel('Velocity ($\mu m/sec$)', fontdict=font)
ax2.plot(time,v, 'g', label='Change in velocity')
ax2.axhline(np.mean(v), linestyle='--', color='g', label='Mean velocity')
ax2.legend(loc='upper right')
ax2.tick_params('y', colors='g')

plt.title('Binding motors [R = 0.2]', fontdict=font)
plt.grid(); plt.show()

Dx_tip = np.diff(bmd03['x']); Dy_tip = np.diff(bmd03['y'])
DD=np.sqrt((Dx_tip**2)+(Dy_tip**2))
v=DD/(2*dt); Av_vel = np.mean(v)
vSD=np.sum(((v-Av_vel)**2)/(np.size(v)-1)); vSD=np.sqrt(vSD)
time = np.linspace(0,300,150)

fig, ax1 = plt.subplots(dpi=500)
#plt.figure(figsize=(14,10))
ax1.set_xlabel('Time', fontdict=font)
ax1.set_ylabel('Number of binding motors', fontdict=font)
ax1.plot(time,bm03, label='Change in binding motor')
#ax1.axhline(float(bm01.mean()), linestyle='--')
ax1.legend(loc='upper left')

ax2 = ax1.twinx()
ax2.set_xlabel('Time', fontdict=font)
ax2.set_ylabel('Velocity ($\mu m/sec$)', fontdict=font)
ax2.plot(time,v, 'g', label='Change in velocity')
ax2.axhline(np.mean(v), linestyle='--', color='g', label='Mean velocity')
ax2.legend(loc='upper right')
ax2.tick_params('y', colors='g')

plt.title('Binding motors [R = 0.3]', fontdict=font)
plt.grid(); plt.show()

Dx_tip = np.diff(bmd04['x']); Dy_tip = np.diff(bmd04['y'])
DD=np.sqrt((Dx_tip**2)+(Dy_tip**2))
v=DD/(2*dt); Av_vel = np.mean(v)
vSD=np.sum(((v-Av_vel)**2)/(np.size(v)-1)); vSD=np.sqrt(vSD)
time = np.linspace(0,300,150)

fig, ax1 = plt.subplots(dpi=500)
#plt.figure(figsize=(14,10))
ax1.set_xlabel('Time', fontdict=font)
ax1.set_ylabel('Number of binding motors', fontdict=font)
ax1.plot(time,bm04, label='Change in binding motor')
#ax1.axhline(float(bm01.mean()), linestyle='--')
ax1.legend(loc='upper left')

ax2 = ax1.twinx()
ax2.set_xlabel('Time', fontdict=font)
ax2.set_ylabel('Velocity ($\mu m/sec$)', fontdict=font)
ax2.plot(time,v, 'g', label='Change in velocity')
ax2.axhline(np.mean(v), linestyle='--', color='g', label='Mean velocity')
ax2.legend(loc='upper right')
ax2.tick_params('y', colors='g')

plt.title('Binding motors [R = 0.4]', fontdict=font)
plt.grid(); plt.show()

Dx_tip = np.diff(bmd05['x']); Dy_tip = np.diff(bmd05['y'])
DD=np.sqrt((Dx_tip**2)+(Dy_tip**2))
v=DD/(2*dt); Av_vel = np.mean(v)
vSD=np.sum(((v-Av_vel)**2)/(np.size(v)-1)); vSD=np.sqrt(vSD)
time = np.linspace(0,300,150)

fig, ax1 = plt.subplots(dpi=500)
ax1.set_xlabel('Time', fontdict=font)
ax1.set_ylabel('Number of binding motors', fontdict=font)
ax1.plot(time,bm05, label='Change in binding motor')
ax1.legend(loc='upper left')

ax2 = ax1.twinx()
ax2.set_xlabel('Time', fontdict=font)
ax2.set_ylabel('Velocity ($\mu m/sec$)', fontdict=font)
ax2.plot(time,v, 'g', label='Change in velocity')
ax2.axhline(np.mean(v), linestyle='--', color='g', label='Mean velocity')
ax2.legend(loc='upper right')
ax2.tick_params('y', colors='g')

plt.title('Binding motors [R = 0.5]', fontdict=font)
plt.grid(); plt.show()

Dx_tip = np.diff(bmd06['x']); Dy_tip = np.diff(bmd06['y'])
DD=np.sqrt((Dx_tip**2)+(Dy_tip**2))
v=DD/(2*dt); Av_vel = np.mean(v)
vSD=np.sum(((v-Av_vel)**2)/(np.size(v)-1)); vSD=np.sqrt(vSD)
time = np.linspace(0,300,150)

fig, ax1 = plt.subplots(dpi=500)
ax1.set_xlabel('Time', fontdict=font)
ax1.set_ylabel('Number of binding motors', fontdict=font)
ax1.plot(time,bm06, label='Change in binding motor')
ax1.legend(loc='upper left')

ax2 = ax1.twinx()
ax2.set_xlabel('Time', fontdict=font)
ax2.set_ylabel('Velocity ($\mu m/sec$)', fontdict=font)
ax2.plot(time,v, 'g', label='Change in velocity')
ax2.axhline(np.mean(v), linestyle='--', color='g', label='Mean velocity')
ax2.legend(loc='upper right')
ax2.tick_params('y', colors='g')

plt.title('Binding motors [R = 0.6]', fontdict=font)
plt.grid(); plt.show()

Dx_tip = np.diff(bmd07['x']); Dy_tip = np.diff(bmd07['y'])
DD=np.sqrt((Dx_tip**2)+(Dy_tip**2))
v=DD/(2*dt); Av_vel = np.mean(v)
vSD=np.sum(((v-Av_vel)**2)/(np.size(v)-1)); vSD=np.sqrt(vSD)
time = np.linspace(0,300,150)

fig, ax1 = plt.subplots(dpi=500)
ax1.set_xlabel('Time', fontdict=font)
ax1.set_ylabel('Number of binding motors', fontdict=font)
ax1.plot(time,bm07, label='Change in binding motor')
ax1.legend(loc='upper left')

ax2 = ax1.twinx()
ax2.set_xlabel('Time', fontdict=font)
ax2.set_ylabel('Velocity ($\mu m/sec$)', fontdict=font)
ax2.plot(time,v, 'g', label='Change in velocity')
ax2.axhline(np.mean(v), linestyle='--', color='g', label='Mean velocity')
ax2.legend(loc='upper right')
ax2.tick_params('y', colors='g')

plt.title('Binding motors [R = 0.7]', fontdict=font)
plt.grid(); plt.show()

Dx_tip = np.diff(bmd08['x']); Dy_tip = np.diff(bmd08['y'])
DD=np.sqrt((Dx_tip**2)+(Dy_tip**2))
v=DD/(2*dt); Av_vel = np.mean(v)
vSD=np.sum(((v-Av_vel)**2)/(np.size(v)-1)); vSD=np.sqrt(vSD)
time = np.linspace(0,300,150)

fig, ax1 = plt.subplots(dpi=500)
ax1.set_xlabel('Time', fontdict=font)
ax1.set_ylabel('Number of binding motors', fontdict=font)
ax1.plot(time,bm08, label='Change in binding motor')
ax1.legend(loc='upper left')

ax2 = ax1.twinx()
ax2.set_xlabel('Time', fontdict=font)
ax2.set_ylabel('Velocity ($\mu m/sec$)', fontdict=font)
ax2.plot(time,v, 'g', label='Change in velocity')
ax2.axhline(np.mean(v), linestyle='--', color='g', label='Mean velocity')
ax2.legend(loc='upper right')
ax2.tick_params('y', colors='g')

plt.title('Binding motors [R = 0.8]', fontdict=font)
plt.grid(); plt.show()

Dx_tip = np.diff(bmd09['x']); Dy_tip = np.diff(bmd09['y'])
DD=np.sqrt((Dx_tip**2)+(Dy_tip**2))
v=DD/(2*dt); Av_vel = np.mean(v)
vSD=np.sum(((v-Av_vel)**2)/(np.size(v)-1)); vSD=np.sqrt(vSD)
time = np.linspace(0,300,150)

fig, ax1 = plt.subplots(dpi=500)
ax1.set_xlabel('Time', fontdict=font)
ax1.set_ylabel('Number of binding motors', fontdict=font)
ax1.plot(time,bm09, label='Change in binding motor')
ax1.legend(loc='upper left')

ax2 = ax1.twinx()
ax2.set_xlabel('Time', fontdict=font)
ax2.set_ylabel('Velocity ($\mu m/sec$)', fontdict=font)
ax2.plot(time,v, 'g', label='Change in velocity')
ax2.axhline(np.mean(v), linestyle='--', color='g', label='Mean velocity')
ax2.legend(loc='upper right')
ax2.tick_params('y', colors='g')

plt.title('Binding motors [R = 0.9]', fontdict=font)
plt.grid(); plt.show()

Dx_tip = np.diff(bmd10['x']); Dy_tip = np.diff(bmd10['y'])
DD=np.sqrt((Dx_tip**2)+(Dy_tip**2))
v=DD/(2*dt); Av_vel = np.mean(v)
vSD=np.sum(((v-Av_vel)**2)/(np.size(v)-1)); vSD=np.sqrt(vSD)
time = np.linspace(0,300,150)

fig, ax1 = plt.subplots(dpi=500)
ax1.set_xlabel('Time', fontdict=font)
ax1.set_ylabel('Number of binding motors', fontdict=font)
ax1.plot(time,bm10, label='Change in binding motor')
ax1.legend(loc='upper left')

ax2 = ax1.twinx()
ax2.set_xlabel('Time', fontdict=font)
ax2.set_ylabel('Velocity ($\mu m/sec$)', fontdict=font)
ax2.plot(time,v, 'g', label='Change in velocity')
ax2.axhline(np.mean(v), linestyle='--', color='g', label='Mean velocity')
ax2.legend(loc='upper right')
ax2.tick_params('y', colors='g')

plt.title('Binding motors [R = 1.0]', fontdict=font)
plt.grid(); plt.show()